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Q ■ Abstract 

^: 

Multiple quantum coherences are typically characterised by their coherence num- 

c , 

^ , ber and the number of spins that make up the state, though only the coherence 

number is normally measured. We present a simple set of measurements that ex- 
^ , tend our knowledge of the multiple quantum state by recording the coherences in 

' both the x basis and the usual z basis. The coherences in the two bases are related 

by a similarity transformation. We characterize the growth of the multiple quan- 
tum coherences via measurements in the two bases, and show that the rate varies 
with the coefficient of the driving term in the Hamiltonian. Such measurements in 
non-commuting bases provides additional information over the ID method about 
the state of the spin system. In particular the measurement of coherences in a basis 
other than the usual z basis allows us to study the dynamics of the spin system 
under Hamiltonians, such as the secular dipolar Hamiltonian, that conserve z basis 
coherence number. 



Preprint submitted to Chemical Physics Letters 



1 February 2008 



Introduction 



The many-body behaviour of nuclear spins in a sohd was described, for a long 
time, in the language of spin thermodynamics [1,2,3]. This description is, how- 
ever, essentially static, and the dynamical behaviour of the system was typi- 
cally ignored, or addressed in terms of memory functions [4,5]. With the advent 
of multiple quantum NMR techniques [6,7,8,9], multi-spin processes could be 
described by multiple spin correlations and multiple quantum (MQ) coher- 
ences. The selective excitation and transformation [10] of MQ coherences led 
to a new picture of many-body spin dynamics in a dipolar solid [11,12,13,14]. 
The focus in these experiments was on transitions of coherence number, which 
were observable, rather than the number of spins involved in the MQ coherence 
states. However, the object of the experiment continued to be "spin counting" 
as these MQ experiments were called. These techniques have been used to 
probe the spatial relationships between spins in large macromolecules, poly- 
mers and crystalline systems, including determining the dimensionality and 
size of localized spin clusters (see [15,16,17,18] for reviews). 

In principle, all the spins in a solid are coupled together through their dipolar 
fields, and the Hilbert space of the spin system is determined by the total 
number of spins in the sample. However, at high field and for temperatures 
above a few degrees Kelvin, it is sufficient to consider a much smaller spin 
system to predict the NMR spectrum of a sohd [19]. The sample resembles 
an ensemble of weakly coupled subsystems, within each of which an effective 
number of coupled spins is postulated to exist. In equilibrium, at high field, 
the spin number is one, as p ~ J^. 
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In a strong magnetic field (Bqz), an N spin-1/2 system has 2^ stationary 
states. These can be classified according to the magnetic quantum number, 
Mz — J2i'>TT'zi — (^|+i/2> — ?^|-i/2>)/2, where m^j = ±1/2 is the eigenvalue of 
the ith spin in the system, and the energy eigenvalue corresponding to is 
Ez = —'jkBoMz- For non-degenerate stationary states there are on the order 
of 2^^~^ possible transitions between any two levels. The difference in 
values between the two levels is referred to as the coherence number. If the 
density operator is expanded in the basis of irreducible tensor operators T;^, 

P = ^aimTim (1) 

Im 

the rank / of the tensor element defines the spin number, while the order m 
characterizes the coherence number. 

While these coherences refer to transitions between levels, it is useful to dis- 
cuss multiple quantum coherences for states of a system. When the state is 
expressed in the eigenbasis of the system, the presence of a non-zero matrix 
element < Zi\p\zj >, indicates the presence of an n quantum coherence, where 
n = Mz{zj) — Mz{zi). Since = J2i''TT'zi is a good quantum number, we use 
a collective rotation about the axis of quantization, J2i to characterize it 

< Zi\exp{—i(l)^Il)pexp{+i(l)^Il)\zj >= exp(m0) < Zi\p\zj > . (2) 

i i 

In the usual MQNMR experiment, the system (initially in Zeeman equilib- 
rium) is allowed to evolve under the action of a Hamiltonian that generates 
single quantum (SQ) [20], or double quantum (DQ) [10,21] transitions. This 
progressively increases the coherence numbers of the state of the system, as 
well as causing its spin number to increase. 
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While the coherences have a physical meaning in the eigenbasis of the sys- 
tem, a generalized coherence number reports on the response of the system 
to any collective rotation of the spins (about the x axis for example). This 
is equivalent to expressing the state of the spins in a basis where the appar- 
ent axis of quantization is given by the axis of rotation, and can be obtained 
from the eigenbasis via a similarity transformation. For example, the similar- 
ity transform V connects the density matrices of the system in the two (the z 
or eigen-basis, and the x basis) representations. 

[p1 = [p^] V (3) 



where the elements of the matrices are =< Xi\p\xj >, [p^]^ =< Zi\p\zj > 
and {xi} and {zi} are complete sets of basis operators. Under a collective 
rotation about the x axis, we obtain, 

< Xi\ exp(— i0^/*)pexp(-|-i0^/*)|a;j >= exp(m0) < Xi\p\xj > (4) 

i i 

where n is the x basis coherence number. Similarity transformations do not 
change the eigen-energies or the physics of the system [22] . Suter and Pearson 
[23] previously used a combination of phase shifts and a variable flip angle pulse 
to encode for coherences in the y basis as well as the z basis. Their technique 
was recently used to study the dynamics of polarization and coherence echoes 
[24] . Requantization in an alternative basis has also been applied to analyzing 
RF gradient NMR spectroscopy [25]. 

In this letter we demonstrate an improved technique for the encoding of coher- 
ences in the x basis as well as for encoding coherences simultaneously in the 
X and z bases. While the measurement of coherence number in an orthogonal 
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basis does provide more information about the state, it does not yield a direct 
measure of the spin number. Under a collective rotation of the spins, the dif- 
ferent orders within a given rank are mixed, but there is no mixing between 
terms of different rank. Thus contributions to a given coherence order from 
the different ranked tensors cannot be separated out, without some measure 
of the distribution of tensor ranks in the system. In order to unambiguosly 
determine the spin number, N independent measurements are required. 

Measurements in non-commuting basis are central to the task of quantum 
state tomography. Eigenbasis measurements provide information on the am- 
plitudes of the terms in the density matrix (in the eigenbasis) , but not on the 
associated phase factors. Changing the basis and repeating the measurements 
allows reconstruction of the exact state of the system, a familiar process used 
to measure the Wigner function in optics experiments [26]. Measuring multi- 
ple quantum coherences in a basis other than the usual z basis is particularly 
important if we wish to study the dynamics of the spin system under a Hamil- 
tonian that conserves z basis coherence number, such as the secular dipolar 
Hamiltonian. 



Methods 



Table 1 shows the initial state, Hamiltonian, and selection rules for the stan- 
dard MQ experiment (using a DQ Hamiltonian) , in both the standard z basis 
and the x basis. Reference [25] tabulates the transformations between quanti- 
zation in the different Cartesian bases. Thus, starting from an initial Zeeman 
state, we see that under the DQ Hamiltonian we should get only even order 
coherences in the z basis and only odd order coherences in the x basis. 
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Table 1 

Description of the MQ experiment in the z and x bases. 





z basis 


X basis 


initial state 


h 


_l //+ _ /-\ 
2 l a; X 5 


initial coherence number 





±1 


MQ Hamiltonian 




i<j 

2 L ' XI XJJ \ 


coherence number selection rule 


±2 


0,±2 


spin number selection rule 


±1 


±1 


coherences 


even 


odd 



The experimental methods presented here improve on those of Suter and Pear- 
son, as their variable flip angle pulse is replaced by a sequence of phase shifted 
pulses, whose duration is fixed. The dipolar evolution during the variable an- 
gle pulse, as it is sampled out to multiples of 27r, can significantly attenuate 
the signal and compromise the resolution of the coherences in the x or y basis, 
especially in a strongly dipolar coupled system. In our experiment, the dipolar 
evolution is refocused, and the 04 rotation achieved purely with phase shifts. 

The pulse sequences shown in Figures 1(a) and 1(b) allow us to encode co- 
herences in the two bases under essentially identical conditions. Figure 1(a) 
is a 2; basis encoding experiment while Figure 1(b) is an x basis encoding 
experiment. The only difference between the two is that the first 7i/2 pulse is 
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phase shifted along with C/^ in the x basis experiment. Figure 1(c) shows the 
16 pulse DQ selective sequence used. It consists of two cycles of the standard 
8 pulse sequence, phase shifted by tt with respect to each other. The sequence 
compensates for pulse imperfections and resonance offsets. C/<^ is created by 
phase shifting all the pulses in the 16 pulse experiment by 0. The two (7r/2) 
pulses and the Cory 48-pulse sequence are not required for the z basis exper- 
iment. However, they are included in order to perform the two experiments 
under identical conditions. In the x basis experiment, the two (7r/2) pulses 
perform the basis transformation, and the phase encoding of the coherences. 
Placing them back-to-back leads to unwanted switching transients, so they are 
separated by the Cory 48-pulse sequence which prevents evolution of the spin 
system under the secular dipolar Hamiltonian between the two (7r/2) pulses, 
and has been described previously [27]. 

The operator corresponding to the observable signal is I^- The measured signal 
for the experiment in Figure 1(a) corresponds to < >(t>— Tr[pfl^, where 
the final density matrix is given by 



where Ra{4>) — exp(i0/Q,), and we have used the fact that the initial state 
is invariant to 2;-rotations. Defining ps — UpiU^ — UIzU\ the state of the 
system after evolution under the DQ Hamiltonian, we obtain the measured 
signal in the z basis experiment 



Pf = C/ti?,(-7r/2)i?,(7r/2)t/<^p,t/ji?,(-7r/2)i?,(7r/2)t/ 



U^R,(-(l>)UhU^RMU 



(5) 



< >0= Tr [Rz{-(f))psRzi(f))ps] 



(6) 
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f X XX XX XX XXXX XX XX x\ 

^|a |a|a |a|a |a|a |a|a |a|a |a|a |a|a | A 

\ A' = 2A + tj^/2 / 
* tc=24(A+W2) * 



Fig. 1. (a) The z basis encoding experiment, = Rz{—4>)URz{4'), and 
U = e^p{iHDQr). The propagator for the 48 pulse time-suspension sequence is 
the Identity Operator I. (b) The x basis encoding experiment where the first (7r/2) 
pulse is phase shifted, (vr/2)(^ = Rz{—(p){TT/2)yRz{4i). (c) The 16 pulse sequence used 
to generate the effective DQ Hamiltonian; A = 1.3 /is, = 0.51 /xs, tc = 43.4 /is. 

For the experiment in Figure 1(b), the final density matrix is given by 

Pf = U^Ry{-7r/2)R^{7r/2)U^p,UlR^i-7r/2)Ry{7r/2)U 

^U^Ryi-7i/2)Rz{-(l>)Ry{7r/2)UIzU^Ry{-n/2)Rz{(j>)Ry{7r/2)U 
^U^R,{-<P)UIzWR,{<P)U (7) 
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and the observed magnetization in the x basis experiment is 
< h >0= Tr [R^{-(t))psRx{(t>)Ps] ■ 



(8) 



In both cases, the experiment is repeated multiple times as is uniformly 
sampled out to a multiple of 27r, and the resulting data Fourier transformed 
with respect to to obtain the distribution of coherence numbers. 

Note that if Equations (6) and (8) could be written in terms of < Iz ><j,= 
Tr[A[(t))ps], where the set of operators A{(f)) form a complete basis for the 
Hilbert space of the spin system, it would be possible to perform quantum 
state tomography on the spins [26]. However, the experiments described here 
cannot completely characterize the state, or even just its collective properties. 

Results 

The experiments were performed at 2.35 T with a Bruker Avance spectrometer 
and a home-built RF probe. The 90 degree pulse time was 0.51 /xs. The pulse 
spacing A used in the DQ sequence was 1.3 /xs, and the cycle time for the 
16 pulse cycle was 43.4 jis. The pulse spacing in the 48 pulse time suspension 
sequence was 1.5 fis. The Ti of the single crystal calcium fluoride sample used 
was 7 s, and the recycle delay used in the experiment was 10 s. 

Figure 2 shows the results obtained in the z and x basis encoding experi- 
ments. The maximum coherence encoded was ±32, with = 27r/64. The 
phase incrementation was carried out to 87r. It is seen that the z and x ba- 
sis measurements give only even and odd coherences respectively as expected 
from Table 1. The data shown correspond to 1, 3 and 5 loops of the 16 pulse 
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Fig. 2. Comparison of z basis and x basis coherences at preparations times r = 
43.4, 130.3 and 217.2 ^us, corresponding to 1, 3 and 5 loops of the 16 pulse cycle 
in Figure 1(c), showing the presence of purely even and odd coherences in the two 
bases respectively. 

cycle. The higher order coherences are seen to grow in both bases, as the 
system evolves under the DQ Hamiltonian. 



In Figure 3 we plot the effective spin number, obtained by a Gaussian fit to 
the coherence number distributions in the x and z bases {N{t) = 2cr^) as a 
function of r [11]. The fits were performed on the ID data. The variance of 
the X basis measurements is consistently smaller than that of the z basis mea- 
surements. The points appear to lie on a straight line, and a best linear fit has 



10 



a slope of 0.54, which is very close to the value of 0.5 expected from the ratio 
between the double quantum selective terms in the DQ Hamiltonian expressed 
in the two bases, as shown in Table 1. The linear fit was performed on the 
mean spin numbers and Nz without considering the standard deviations. 
The error in the fit is ncghgible. In the z basis, evolution under the DQ Hamil- 
tonian forces the system to change coherence number, while in the x basis the 
presence of zero quantum terms permits mixing without changing the coher- 
ence number. Thus the growth of the coherence numbers is slowed relative to 
the z basis, leading to a narrower distribution. It should be emphasised that 
the basis change does not change the spin number, only the experimentally 
observable coherences. The change in spin number with basis representation 
demonstrates the limitations of the Gaussian statistical model as an accurate 
predictor of spin number in strongly coupled spin systems. Lacelle has also 
discussed the limitations of the model [17]. 

As first shown by Suter and Pearson [23] , a two dimensional experiment illus- 
trates the correlation between x and z basis coherences. The two-dimensional 
experiment is obtained from the x basis experiment in Figure 1(b) by phase 
cycling the refocussing sequence C/^ by (3 independently of 0. The phases (f) 
and (5 are incremented independently to sample a rectangular grid and a 2D 
Fourier transform is performed to yield the coherences. The measured data in 
a single experiment is 

< h >/3</>= Tr [Rz{-P)R,{-(l>)psRMRz{P)Ps] . (9) 

It is straightforward to show that the order of the x and z phase shifts does 
not matter when both of them are sampled over a 27r range. The {(t>x){l3z) 
experiment is equivalent to the {—I3z){—4'x) experiment. The two dimensional 



11 




180 



Fig. 3. Plot of the effective spin number (-/V(r) = 2(7^) obtained by fitting the co- 
herence distributions obtained in the ID x and z basis measurements to a Gaussian. 
Also shown is the best linear fit to the data, whose slope is 0.54. The slope expected 
from the coefficient of the DQ selective terms in the two bases (see Table 1) is 0.5. 

experiment separates out the different terms that contribute to a particular z 
basis coherence as can be seen in Figure 4. We used r = 130.3 //s, corresponding 
to 3 loops of the MQ cycle. The maximum coherence encoded in each direction 
was ±12, with A0 — 27r/24. The phase incrementation was carried out to Stt 
along each axis, resulting in a 96 x 96 data grid, which was Fourier transformed 
to yield the coherences shown. 

This two dimensional technique can be used to examine the evolution of mul- 
tiple quantum coherences under the dipolar Hamiltonian. Figure 5 shows the 
attenuation of the z basis zero quantum signal as it evolves under the dipolar 
Hamiltonian. Also shown on the figure are the various x basis contributions 
to the z basis zero quantum signal obtained from the 2D data. The decay is 
clearly non-exponential. It can be seen that the different x basis terms at- 
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Fig. 4. Results of the 2D experiment - showing correlations between encoding in 
the X and z bases. The preparation time used was r = 130.3 ;us, corresponding to 
3 loops of the 16 pulse cycle in Figure 1(c). The width in z appears broader than 

the width in x. 



tenuate at different rates, and that the measured decay of the ID z basis 
data represents some sort of average attenuation of all these terms. Thus this 
technique allows us to probe the details of spin dynamics beyond the abil- 
ity of existing techniques. We have also recently used this technique to study 
the evolution of the spin system following a Jeener-Broekaert pulse pair, and 
observed the evolution of the system to a dipolar ordered state [28] . 
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Fig. 5. Decay of the z basis zero quantum signal under the dipolar Hamiltonian. 
The decay of the different x basis contributions to the zero quantum signal in the z 
basis, obtained from the 2D experiment are also shown. The data were normalized 
to the observed intensity at 2 /xs. Note that the data are plotted on a log scale. The 
preparation time used was r = 130.3 fj,s, corresponding to 3 loops of the 16 pulse 
cycle in Figure 1(c). 



Conclusions 



We have shown that by encoding MQ coherences in different bases {x and z) 
additional information about the state of the spin system may be obtained. 
In particular, x basis encoding could be useful in determining the size of 
multiple spin correlations under the action of a Hamiltonian that preserves z 
basis coherence number, but changes the number of spins in the state, such 
as the secular dipolar Hamiltonian. 
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Figure Captions 

(1) (a) The z basis encoding experiment, — Rz{—(^)URz{<^^-, and U — 
exp{iHDQT). The propagator for the 48 pulse time-suspension sequence 
is the Identity Operator I. (b) The x basis encoding experiment, where 
the first (7r/2) pulse is phase shifted, — Rz{—(l)){n/2)yRz{(l)). (c) 
The 16 pulse sequence used to generate the effective DQ Hamiltonian; 
A = 1.3 fjLS, = 0.51 /is, tc — 43.4 //s. 

(2) Comparison of z basis and x basis coherences at preparations times r = 
43.4, 130.3 and 217.2 /xs, corresponding to 1, 3 and 5 loops of the 16 
pulse cycle in Figure 1(c), showing the presence of purely even and odd 
coherences in the two bases respectively. 

(3) Plot of the effective spin number {N{t) — 2a^) obtained by fitting the 
coherence distributions obtained in the ID x and z basis measurements 
to a Gaussian. Also shown is the best linear fit to the data, whose slope 
is 0.54. The slope expected from the coefficient of the DQ selective terms 
in the two bases (see Table 1) is 0.5. 

(4) Results of the 2D experiment - showing correlations between encoding 
in the x and z bases. The preparation time used was r = 130.3 /is, 
corresponding to 3 loops of the 16 pulse cycle in Figure 1(c). 

(5) Decay of the z basis zero quantum signal under the dipolar Hamiltonian. 
The decay of the different x basis contributions to the zero quantum 
signal in the z basis, obtained from the 2D experiment arc also shown. 
The data were normalized to the observed intensity at 2 /is. Note that 
the data are plotted on a log scale. The preparation time used was r = 
130.3 /is, corresponding to 3 loops of the 16 pulse cycle in Figure 1(c). 
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